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Abstract 

We compute the helicity skewed quark distributions H and E in the chiral quark-soliton 
model of the nucleon. This model emphasizes correctly the role of spontaneously broken 
chiral symmetry in structure of nucleon. It is based on the large- A^c picture of the nucleon 
as a soliton of the effective chiral lagrangian and allows to calculate the leading twist quark- 
and antiquark distributions at a low normalization point. 

We discuss the role of chiral symmetry in the helicity skewed quark distributions H and 
E. We show that generalization of soft pion theorems, based on chiral Ward identities, leads 
in the region of — ^ < a: < ^ to the pion pole contribution to E which dominates at small 
momentum transfer. 
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1 Introduction 



Recently, a new type of parton distributions [||, ||, ^ ^ has attracted considerable interest, 
the so-called skewed parton distributions (SPD's), which are generalizations simultaneously of 
the usual parton distributions, distribution amplitudes and of the elastic nucleon form factors 
(for review see [^). Taking the n-th moment of the SPD's one obtains the form factors (i.e., 
non-forward matrix elements) of the spin-n, twist-two quark and gluon operators. On the 
other hand, in the forward limit the SPD's reduce to the usual quark, antiquark and gluon 
distributions. In other words, the SPD's interpolate between the traditional inclusive (parton 
distributions) and exclusive (form factors) characteristics of the nucleon and thus provide us 
with a considerable new amount of information on nucleon structure. 

The SPD's are not accessible in standard inclusive measurements. They can, however, be 
measured in deeply-virtual Compton scattering (DVCS) and in hard exclusive leptoproduction 
of mesons. The very possibility to probe SPD's in these reactions is due to QCD factoriza- 
tion theorem of Ref. Q. Feasibility of experimental measurements of SPD's in hard exclusive 
reactions is currently being studied 0, H, |9|. A quantitative description of these classes of 
processes requires non-perturbative information in the form of the SPD's at some initial nor- 
malization point. Although the skewed parton distributions can be reduced in certain limiting 
cases to already known quantities (parton distributions, form factors), even their qualitative 
behaviour is unknown to large extent. That is why model calculations of these quantities are 
of big importance. There were already model calculations of SPD's: in the bag model and 
in the chiral quark-soliton model In the latter calculation a drastic variation of flavour 

singlet H{x,^,t) at x near x = |^| was observed. Such behaviour is related to the fact that the 
SPD's in the region — ^ < a; < ^ have properties of distribution amplitudes. This feature, being 
very important for the understanding of SPD's, requires a field theoretic description of nucleon's 
constituents and that is the reason why it can not by reproduced in the bag model. 

Our aim now is to compute helicity skewed quark distributions of the nucleon using the 
methods of ref. |11]. We shall see that generalization of low energy theorems requires that the 
skewed distribution E develops a pion poleQ at = of the form: 



where $7r(-z) is distribution amplitude of the pion. In refs. [^, |l^] it was shown that this contri- 
bution to E leads to considerable enhancement of the amplitude of hard exclusive production 
of charged pions and to large azimuthal spin asymmetry in exclusive tt^ production p^ . 

We shall see that in the chiral quark-soliton model the pion pole contribution is related to 
the large distance asymptotic of the pion mean-field, which is controlled by PCAC. 



2 Definition of skewed helicity quark distributions 

In QCD the helicity skewed quark distributions are defined through non-diagonal matrix ele- 
ments of product of quark fields at light -cone separation. Here and in the following, we shall 
use the notations of ref. ||5[ 

^ The contribution of the pion pole to E was discussed at quaUtaive level in ref. |12] 
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+ ^Ej{x,C,A^)U{P'){n.Ah,U{P). 



(2.2) 

Here is a light-cone vector, 

= 0, n-{P + P') = 2n-P = 2, (2.3) 

A is the four-momentum transfer, 

A = P' -P, (2.4) 

Mtv denotes the nucleon mass, and U {P) is a standard Dirac spinor. The skewed quark distri- 
butions, i/(x,^, A^) and ^^(x,^, A^), are regarded as functions of the variable x, the square of 
the four-momentum transfer, A^ = t, and its longitudinal component 

e = -^(n.A). (2.5) 

In the forward case, P = P' , both A and ^ are zero, and the second term on the r.h.s. of 
eq.( |2.2D disappears. In this limit the function H becomes the usual polarised parton distribution 
function. 

On the other hand, taking the first moment of eq.( |2.2[ ) one reduces the operator on the l.h.s. to 
the local axial vector current. The dependence of H and E on disappears, and the functions 
reduce to the usual axial form factors of the nucleon, 

'dxH{x,tA^) = Ga{A^), (2.7) 

-1 
1 _ 

dxE{x,C,A'^) = Gp(A2). (2.8) 

1 

Taking higher moments of the distribution functions one obtains the form factors of the 
twist-2, spin-n operators. 



3 Chiral quark— soliton model of the nucleon 



Recently a new approach to the calculation of quark distribution functions of the nucleon has 



been developed |14] in the framework of the chiral quark-soliton model of the nucleon [15|. In 
present paper we apply this approach to the calculation of skewed quark distributions. It is 
essentially based on the 1/Nc expansion. Although in reality the number of colours Nc = 3, the 
academic limit of large is known to be a useful guideline. At large the nucleon is heavy 
and can be viewed as a classical soliton of the pion field |16, 17|. In this paper we work with the 
effective chiral action given by the functional integral over quarks in the background pion field 
MM: 
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exp(iS'efr[^(a;)]) = / D^jDiIj exp (i f (fxi^iid - MW-)^) 



U = exp(i7r'^(x)T"), ?7^5 = exp (i^'^(x)r"75) = + (3.9) 

Here is the quark field, M is the effective quark mass, which is due to the spontaneous break- 
down of chiral symmetry (generally speaking, it is momentum dependent), and U is the SU{2) 
chiral pion field. The effective chiral action given by ( |3.9| ) is known to contain automatically the 
Wess-Zumino term and the four-derivative Gasser-Leutwyler terms, with correct coefficients. 
The equation ( |3.9| ) has been derived from the instanton model of the QCD vacuum |2C, ^ ], 



which provides a natural mechanism of chiral symmetry breaking and enables one to express 
the dynamical mass M and the ultraviolet cutoff intrinsic in ( |3.9| ) through the Aqqd parameter. 
The ultraviolet regularization of the effective theory is provided by the specific momentum de- 
pendence of the mass, M(p^), which drops to zero for momenta of order of the inverse instanton 
size in the instanton vacuum, 1/p ~ 600 MeV. For simplicity we shall neglect this momen- 
tum dependence in the general discussion; it will be taken into account again in the theoretical 
analysis and in the numerical estimates later. 

An immediate application of the effective chiral theory (|3.9| ) is the quark-soliton model 
of baryons of ref. |15], which is in the spirit of the earlier works [ p2| , p^ . According to this 



model nucleons can be viewed as Nc "valence" quarks bound by a self-consistent pion field (the 
"soliton" ) whose energy coincides with the aggregate energy of the quarks of the negative-energy 
Dirac continuum. Similarly to the Skyrme model large Nc is needed as a parameter to justify 
the use of the mean-field approximation; however, the 1/A^c-corrections can be — and, in some 
cases, have been — computed j2j]. 

Let us remind the reader how the nucleon is described in the effective low-energy theory 
(|3.9|). Integrating out the quarks in (|3.9|) one finds the effective chiral action. 



5effK(x)] = -iVeSplogZ)(C/), D{U) = ido-H{U), (3.10) 

where H(U) is the one-particle Dirac Hamiltonian, 

H{U) = -i-f'^-f^dk + Mj^U^^ , (3.11) 

and Sp . . . denotes the functional trace. For a given time-independent pion field U = exp(i7r"(x)r") 
one can determine the spectrum of the Dirac Hamiltonian, 

H<^>n = EA. (3.12) 

It contains the upper and lower Dirac continua (distorted by the presence of the external pion 
field), and, in principle, also discrete bound-state level(s), if the pion field is strong enough. If the 
pion field has unity winding number, there is exactly one bound-state level which travels all the 
way from the upper to the lower Dirac continuum as one increases the spatial size of the pion field 
from zero to infinity |l^. We denote the energy of the discrete level as E\ev, —M < Ei^y < M. 
One has to occupy this level to get a non-zero baryon number state. Since the pion field is colour 
blind, one can put Nc quarks on that level in the antisymmetric state in colour. 

The limit of large Nc allows us to use the mean-field approximation to find the nucleon 
mass. To get the nucleon mass one has to add A^c-f'iev ^-nd the energy of the pion field. Since 
the effective chiral lagrangian is given by the determinant ( |3.10| ) the energy of the pion field 
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coincides exactly with the aggregate energy of the lower Dirac continuum, the free continuum 
subtracted. The self-consistent pion field is thus found from the minimisation of the functional 
11 

Mn = min N, \ Ei,^[U] + J2 {^n[U] - 4°^) [ • (3-13) 

From symmetry considerations one looks for the minimum in a hedgehog ansatz: 

f/c(x) = exp (i7r"(x)r'') = exp (m''r''P(r)) , r = |x|, " = ~' (^"^^^ 

where P{r) is called the profile of the soliton. 

The minimum of the energy ( 3.13| ) is degenerate with respect to translations of the soliton 



in space and to rotations of the soliton field in ordinary and isospin space. For the hedgehog 
field ( |3.14| ) the two rotations are equivalent. The projection on a nucleon state with given spin 
(53) and isospin (Ts) components is obtained by integrating over all spin-isospin rotations, R 

{S = T, S3, n\...\S = T, S3, Ta) = JdR (p^^jfiR) ■ ■ ■ 47n{R) ■ (3-15) 

Here (f'sTT-A-^) is the rotational wave function of the nucleon given by the Wigner finite-rotation 
matrix [0, 

CtI(^) = i-lf^'^^V2STTD%%{R). (3.16) 

Analogously, the projection on a nucleon state with given momentum P is obtained by integrat- 
ing over all shifts, X, of the soliton, 

(P'|...|P) = y d^X e^(^'-P)-^ . . . (3.17) 

4 Skewed quark distributions in the chiral quark— soliton model 

We now turn to the calculation of the skewed quark distributions in the chiral quark-soliton 
model. This description of the nucleon is based on the l/A'^c-expansion. At large Nc the nucleon 
is heavy — its mass is 0{Nc). For the large- A'c nucleon eq. ( ^.2| ) simplifies as follows: 



^e^^^P',S!,\iJf{-Xn/2)h^r,^l^fiXn/2)\P,S3) 



2 ^25'^Hfix,C,t) - -^Ep{x,^,t)^a^siS,^ (4-18) 



where 5*3, 5*3 denote the projections of the nucleon spin. From this expression we immediately 
see that in the leading order of the l/A'c-expansion only the flavour isovector part of 

H^^\x, A^) = H^{x, A2) - Hd{x, A2) 

and 



(a;, a2) = E^{x, A2) - Edix, A 



2^ 
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are non-zero. The isosinglet part of H{x,^, A'^) and the isosinglet part of E{x,^, A'^) appear 
only in the next -to-leading order of the l/A'c-expansion, i.e., after taking into account the finite 
angular velocity of the soliton rotation. 

Before computing the skewed quark distribution functions we must determine the parametric 
order in l/N^ of the kinematical variables involved. Generally, when describing parton distri- 
butions in the large- A'^c limit, one has x ~ 1/Nc, since the nucleon momentum is distributed 
among Nc quarks. Furthermore, as in the calculation of nucleon form factors we consider mo- 
mentum transfers to be of order t ~ N^; hence, in particular, ^ ~ l/-^c) so that ^ is of the same 
parametric order as x. 

Technically the calculation of the skewed parton distributions proceeds in much the same 
way as that of the usual parton distributions [0, |ll| . Using the formalism developed in |ll| 
we obtain: 



m^\x,t^^) = - 1 dz" 1 d^X exp (i A • X) J2 exp(izO[(x + OM;v-i?n 

^ occup. 

cDt (X)(AiT3 + 2^MnA^ ■ Tx) (l + 7V)75 (X - zOea), (4.19) 



E^^\x,^,A^ 



NrMl 



dz^ j (fx exp (iA • X) exp (iz^ [{x + ^Mn - En 



?,TTiA\ 

4(X) ( ■ ) ( 1 + 7V )75 (X - z^es). 



occup. 



(4.20) 



Before going ahead with the evaluation of the expressions eqs.( [4.19| , 4.20) we would like to 
demonstrate that the two limiting cases of the skewed distributions — usual parton distributions 
and elastic form factors — are correctly reproduced within the chiral quark-soliton model. 



Taking in eg. ( 4. 19 ) the forward limit, A — > 0, one recovers the formula for the usual polarized 



(anti-) quark distributions in our model which was obtained in ref. [14|. Thus the forward limit, 
eq. (p.6|), is reproduced. On the other hand, integrating eq. (4.19) over — 1 < x < 1 one obtains 
(up to corrections parametrically small in 1 /Nc) the expressions for the axial form factors of the 
nucleon derived in ref. |25|: 



j\xH^^\x,i,A^) = -^|d3Xexp(iA-x) ^ ci>1^(X) rV7S (X) 



occup. 



G' 



(T=l) 



(A2 



(4.21) 



Actually experimental is very well reproduced in the chiral quark-soliton model up 



to momenta of order ~ 1 GeV^ |2£^ 

Now if one integrates eq. ( f4.20| ) over — 1 < x < 1 one obtains (up to corrections parametrically 
small in 'i/Nc) the following expression: 



dxE^^\x,i,A^) 



2NcMn 



J d^Xexp (iA • X) $]^(X)7^'(A^ • Ti)75<&n(X). 



occup. 



(4.22) 

Using the "hedgehog" symmetry of the pion mean-field one can easily show that the expression 
( 4.22 ) is a function of only A^ and coincides with expression for pseudoscalar nucleon form factor 



in the chiral quark soliton model, see e.g. ref. |24]. 
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Eqs. ( [4.19 , 4.20| ) express the SPD's as a sum over quark single-particle levels in the soliton 
field. This sum runs over all occupied levels, including both the discrete bound-state level and 
the negative Dirac continuum. We remind the reader that in the case of usual parton distribu- 
tions it was demonstrated that in order to ensure the positivity of the antiquark distributions it 
is essential to take into account the contributions of all occupied levels of the Dirac Hamiltonian 
1 14]. We shall see below that also in the case of skewed quark distributions the contribution of 
the Dirac continuum drastically changes the shape of the distribution function. That is espe- 
cially important to reproduce the pion pole contribution to the spin-flip SPD E^^^ required by 
chiral Ward identities. 

The contribution of the discrete bound-state level to eqs. ( 4.19 , 4.2[1| ) can be computed using 
the expressions given in the Appendix. The result is shown in Figs. 1 for the forward case and 
Figs. 2 for a non-zero momentum transfer. Being taken by itself this contribution resembles 
qualitatively the shape of SPD's H and E obtained in the bag model ||lOt| . 

To calculate the contribution of the Dirac continuum to eqs. we resort to an 

approximation which proved to be very successful in the computation of usual parton distribu- 



tions, the so-called interpolation formula [14|. One first expresses the continuum contribution 
as a functional trace involving the quark propagator in the background pion field. The quark 
propagator can then be expanded in powers of the formal parameter dU /(— 9^ + M"^), which be- 
comes small in three limiting cases: i) low momenta, \dU\ <C M, ii) high momenta, \dU\ ^ M, 
in) any momenta but small pion fields, |logC/| <C 1. One may therefore expect that this ap- 
proximation has good accuracy also in the general case. As was shown in refs. [14| for usual 
parton distributions this approximation preserves the positivity of the antiquark distributions 
and all sum rules; moreover, it gives results very close to those obtained by exact numerical 
diagonalisation of the Dirac Hamiltonian and summation over the negative-energy levels. 

Simple generalization of the technique developed in ^] allows us to express the Dirac 
continuum contribution to the SPD's directly in terms of pion mean-field ( p. 14 ): 



^cont 



MnNc 
M- 



Im 



{2ttY J {2tt) 



^ Mat — V ■ p 



1/2 



3/2 



(^pj - Mf + io) - M| + iO 



k-v + -A 



(4.23) 



X Trfl. 



+ 



e'-^\x,C,A'^ 



/cont 



2MlN, 



3? Ai 



Im 



d^k f d^p 



(2^)3 7 (2^) 



M- 



1/2 



3/2 



(pj - Mf + io) - M| + iO 



(^x — C^Mat — V ■ p 
[k-v + -A-v^ 



(4.24) 



X Trfl. 



(A^.r^)c/(k-f)[Lr(k+f) 



+ k--e, A^-A , 



where v = (1,0,0,-1) is a light cone vector and the Fourier transform of the soliton field is 
defined as 
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Uik) = j cfxe-'^-''U{x). (4.25) 
Also we introduced short notation pi = p + A, p2 = p — k — A/2, Mi = M{p\) and M2 = M{p2)- 



In eqs. (4.23, 4.24) the momentum dependence of the constituent quark mass, M(p^) cuts 
the loop momentum p and thus regularizes the UV divergence. Let us note that expressions 
( 4.23| , 4.24) are explicitly symmetric under transformation ^ what follows from charge 

conjugation symmetry |5|. 



4.1 Pion pole contribution to the skewed quark distribution E 

Before presenting the numerical results for the SPD's H and E let us discuss specific contribution 
to these SPD's originating from the long range pion tail of the pion mean-field. The behaviour of 
the mean pion field at large distances is governed by linearized equations of motion and PCAC: 



Jim U (x) 



8vr/2|x|2 



(4.26) 



where qa ~ 1-25 is the axial charge of the nucleon, f^^ ~ 93 MeV is the pion decay constant. 
This asymptotic implies that the Fourier transform ( 4.25| ) has the following small momentum 
asymptotic: 



Uik) ~ i2TTf6{k) + 



^9A (k-r) 



+ mi. 

It is useful to split the Fourier transform of the pion mean-field into two pieces: 



(4.27) 



U{k) = C/(A:)s„.ooth + {2TTf 5{k) , 



where 



U{k) 



smooth 



(4.28) 



(4.29) 



Now if we substitute the representation of the Fourier transform of the pion mean- field ( 4.2g| ) 
into expressions ( 4.23| , 4.24) we see immediately that the delta function piece in eq. (4.28) does 
not contribute to H, which means that in (|]2|) we can always replace U by its smooth part 
f^smooth- On contrary in the expression (|4.24| ) for E the contribution of the delta function is 
nonzero and has the form (we denote this contribution Et^): 



Ei^\x,^,A' 



F{A^ 



mV2((p + A)- 



r 



d'^p 
(2^ 

M3/2((P- A) 



(p + A)2-M2-hiO) ((p- A)2 - M2 -KiO 



V ■ p 



+ u 



(4.30) 



where we introduced the following form factor: 

4M2,/2 



F{-k' 



3k^ 



J (fix ex.p{ik ■ x) Tr (u{ 



1 



(4.31) 
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Now the crucial observation is that the integral over p coincides exactly (up to trivial renaming 
of variable) with the expression for the light-cone pion distribution amplitude in the instanton 
model of the QCD vacuum m, Therefore the expression ( |4.30D for the E,^ can be written 
in the compact form: 



where <J>7r(-Z7 A^) is the distribution amplitude of virtual pion, normalized by 



p dz'^^{z,A^) = 1. 



(4.32) 



(4.33) 



Generalizing slightly the technique of refs. |26, 27 1, we computed the dependence of the 
virtual pion distribution amplitude. At small it has the form: 



$^(z,A2) = 1 



ml 



247r2/2 



\^^{z) + ^(1-^') 



N,{A^ - ml) 
247r2/2 



(4.34) 



The pion distribution amplitude calculated in the instanton model of QCD vacuum |26, 27] is 
very close to the asymptotic one ^■k{z) = ^t^{z,A? ml) = |(1 — z^). Therefore from the 
eq. ( ^.34| ) we can conclude that the dependence of the distribution amplitude on the virtuality 
of the pion is rather weak. In addition, under evolution this dependence disappears and we have 
asymptotically '^■niz, A?)asy = |(1 — -s^^)- In what follows we shall therefore drop the dependence 
of distribution amplitude of the virtual pion on A2. 

Using the small momentum asymptotic of the pion mean-field ( 4.27] ) one gets immediately 
the small A2 asymptotic of the form factor F(A2) of eq. ( [4.31 ): 



A2 



lim F(A2 



N 



(A2 



mt 



The expression (4.35) yields then the pion pole contribution to the SPD E'^'^^: 



hm 43)(^^^^^2) ^ 

and as a consequence of the sum rule (| 
form factor G'^p=^ {A'^): 



AgAM_ 



N 



iei(A2 



mt 



■ e(\x\ < |ei)«>. 



X 



(4.35) 



(4.36) 



the pion pole contribution to the pseudoscalar nucleon 



lim Gl=^{A 



A2 



lim 



dxEi'\x,C,A^) 



(A2 -m2 ) 



(4.37) 



We see that the appearance of the pion pole in S^^^ is required by spontaneously broken chiral 
symmetry. In order to reproduce it in some model the latter should respect the chiral Ward 
identities. For example, in computation of SPD's in the bag model pH] the chiral Ward identities 



are violated and the pion pole contribution ( |4.35D is missed. The chiral quark-soliton model 
respects all chiral Ward identities what allows to split unambiguously SPD E^'^'> into two pieces: 



(a;, ^, t) = ^, t) + i^JSooth t) . 



(4.38) 



^More precisely, with distribution amplitude of virtual pion with virtuality 
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the result for E^^ is given by eq. ( 4.32| ), the results for E^'^^^thi^^ given by eq. ( [4.24| ) with 



U replaced by its smooth part [/smooth- 

Let us note that the form factor F{t) (t = A^) in a parametrically wide region m| <C |i| ^ 



contains significant contributions other than the simple pion pole ( 4.35 ). Numerically we 
found that the form factor F{t) can be parametrized at \t\ <C in the following form: 

4.4 / l.l(t-ml)\ 
F{t) w T 1 + ^ , (4.39) 

where all dimensional quantities are in units of GeV. 

It is worth mentioning that although we obtain the pion pole contribution to the spin-flip 
SPD E ( [4. 361 ) in the model calculation, actually the existence of this contribution as such follows 
from general considerations of the chiral Ward identities [0] . 

5 Numerical results and discussion 

We have calculated numerically the isovector distributions ^, A^) and ^, A^). 

For the calculations we use the variational estimate of the soliton profile, eq.( |3.14| ), of ref.[^] 
(Mo = 350 MeV), 

P{r) = -2 arctan , tq ^ I.O/Mq, Mn ~ 1170 MeV, (5.40) 

which has been used in the calculation of usual parton distributions in refs.|14|. Furthermore, 
we approximate the momentum-dependent mass predicted by the instanton model of the QCD 
vacuum ||2^ by the simple form 

^'(V) = (5.41) 

where the parameter A is related to the averaged instanton size, p, by A = 6^/3p-i^ This 
expression reproduces the asymptotic behaviour of Mijp') at large euclidean obtained in the 
instanton vacuum, 

We have explored also other forms of the momentum dependence of the mass and found that 
numerically the results are very close to each other. 

5.1 Results for //(x,^,A2) 

We estimate the Dirac continuum contribution to -ff^^^(x, ^, A^) using the interpolation for- 
mula, eq. ( 4.23| ), which gives a reliable approximation preserving all qualitative features of the 



continuum contribution. The contribution of the discrete level is calculated using eq. ( |A.5| ). 

First we compute H^'^\x,^, Is?) in the forward limit, A ^ 0, where it coincides with the 
usual quark and antiquark distributions: 

#(3)(^,^ = 0,A2 = o) = [ ^^^'^'^VtV'. ^ n' (5-42) 

V ' / I Au(-x) - A(i(-x), X < 0. ^ ' 
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The result is shown in Fig. 1, where we plot separately the contributions of the discrete level and 
that of Dirac continuum (computed from the interpolation formula), as well as their sum. The 
forward limit reproduces the polarized quark distribution obtained in |jl^ in the same model and 
at the same level of approximation. From Fig. 1 we see that the Dirac continuum contribution 
leads to new qualitative prediction for polarized quark distribution: the existence of large flavour 
asymmetry of antiquark distribution Au{x) — Ad{x), the feature which was noted first in |]l^] 
(see also [p^]). To our best knowledge all parametrizations of polarized quark distributions 
assume flavour symmetric antiquark distributions. 

The importance of the Dirac continuum contribution to H^^'^ {x, A^) is also shown on Fig. 2. 
We plot there separately the discrete level and continuum contributions for = —0.5 GeV^ 
and ^ = 0.2. Comparing Fig. 1 and Fig. 2 we see that the ^-dependence of H^^\x,S^,A'^) is 
mostly due to the ^-dependence of Dirac continuum contribution. Since the Dirac continuum 
contribution is symmetric in variable x, we can expect strong ^-dependence only in part of H^^^ 
which is even in x . This observation allows us to conclude that with good accuracy we can put 
(about A^ dependence see below) 

#(3) {x, - H'^^^ {-X, ~ ^u{x) - Au{x) - Ad{x) + Ad{x) . 

The message which might be useful for modelling of SPD's in terms of double distributions 
1^ Additionally we see that "antiquark" skewed distribution {H^^^ at negative x) is large 
and it is dominated by the Dirac continuum contribution. 

In order to illustrate the dependence of H^^\x, ^, A^) on ^ and A-^ we plot this function for a 
fixed momentum transfer of A^ = —0.5 GeV^ for various values of ^ (see upper panel of Fig. 3), 
and for fixed = 0.2 and various values of momentum transfer (see lower Fig. 3). We clearly 
see from the figures that the distribution has "cusps" at x = it^. 

Using the results for H^^\x, ^, A'^) we computed ^ and A-^ dependence of its Mellin moments. 
Lorentz invariance requires that the A^-th Mellin moment should be a polynomial of order in 

lN/2] 

dx x''-'h'^'\x,^,A') = J2 f hki^')- (5.43) 
Here we prefer to use decomposition of Mellin moments in partial waves of qq pairs in t— channel 



i 



dx X 



N-l 



(5.44) 



Where Pi{x) are Legendre polynomials and I is an angular momentum of exchanged qq pair, it 
runs over (odd) even values for (odd) even A^. If we now take the forward limit in eq. ( 5.44| ) we 
obtain: 



dx X 



N-l 



Au{x) - Ad{x) 



1 

2^ 



r(2A^ + 1) 



AN) 



(0). 



(5.45) 



Qualitatively we may expect that the slope of A'^ dependence of the form factor a[^'*(A'^) is 
governed by the mass of a low-lying isovector resonance with spin / and unnatural parity. This 
dependence can be phenomenologically described by simple dipole fit: 
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(Nl) 




M^"^ GeV 


(11) 


1.25 


0.9 


(20) 


0.06 


0.9 


(22) 


0.12 


1.1 


(31) 


0.08 


0.9 


(33) 


0.04 


1.4 



Table 1: Values of parameters of dipole fit for Mellin moments of isovector H, see eq. (5.46) 



AN) 



(0) 



'(A^) = 7 ^^7]^' (5-46) 



where M^^'^ is a phenomenological parameter (dipole mass) . The results of the calculation and of 
the fits are given in Table 1. One should note that a!~^\^) is identical to the usual axial coupling 
constant qa and m}^"* to the corresponding dipole mass parameter. We see that our theoretical 
values of these parameters are very close to experimental ones. One should not overestimate this 
in view of the approximations done in the present approach (interpolation formula, no rotational 
1/A^c corrections, etc.). Still the result shows that the method is well founded. As far as the 
other a\^\^^) concerns, we see that the chiral quark-soliton model reproduces the qualitative 

expectation that the dipole mass in the form factors ai^\l^), describing an exchange with 
angular momentum I, in fact increasing with /. 

In modelling of A^-dependence of SPD's usually the factorization ansatz H{x,S^, /S?) = 
H{x,^)Ga{^^) is used. This ansatz would imply e.g. that dipole masses m/^-* are independent 
of and I. By explicit calculation in the chiral quark-soliton model we demonstrated that this 
is not the case. 

Let us note that the hard exclusive reactions can be viewed as a "tool" to create fundamental 
probes which are absent in the Nature. For example, among the form factors a\^\/!^) only one 
{N = 1, / = 1) can be measured by a probe provided by the Nature {W, Z bosons), all others 
are not accessible for electroweak probes. Since in hard exclusive reactions we can extract the 
higher spin form factors, this allows us to study low energy observables with probes of any spin. 
In this respect the hard exclusive reactions can be used not only for checking of predictions of 
perturbative QCD but also as a new way to study low energy properties of hadrons. 



5.2 Results for E{x, ^, A 



2\ 



The skewed distribution E^^\x,^,A'^) is dominated by the pion pole contribution ( 4.32| ). We 



computed also the smooth part of E^^\x,(,, A"^) (see eq. ( 4.38D ). The results are presented in 



Figure 4, where we plot the pole and the smooth parts separately at various values of ^ and 
A^. We see that the pole contribution dominates E^^^x, A^) in large range of A2 and The 
results for the total distribution (pole+non-pole) at A^ = —0.5 GeV^ and various values of ^ are 
presented on Fig. 5. The ^ dependence of Mellin moments of E^^^ (x, ^, A^) has also polynomial 
form and we checked that our model calculations reproduce this feature. 

The arguments presented here for E can be easily extended to the analogous SPD for N ^ Y 
SPD's (Y is a hyperon from the octet Y = A,T,). In this case the flavour changing E has a 
contribution from of the kaon pole of the form: 
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E- 



.^^y _ 4<-^M^ r X (x 



{\^\<\^\)^k(^], (5.47) 



where is kaon distribution amplitude and the constant entering in description of the 

semileptonic decays of hyperons. 

5.3 Comparison with other models 



Hehcity skewed quark distributions were computed previously in the bag model |10|. Unfortu- 
nately in this model the chiral symmetry is broken explicitly by boundary conditions at bag 
surface. Therefore the crucial contribution of the pion pole to ^, A^) is missed in this 

model. Nevertheless the bag model describes qualitative features of the skewed quark distri- 
butions for which the "resonance part" Q is relatively small, these are odd (even) in x part of 
isovector (isoscalar) H{x,S,,A) and E{x,(,,A). 

In another approach, proposed by Radyushkin ||3|, one writes a spectral representation 
for the matrix element of the light-ray operator in terms of a so-called double distribution. 
The skewed distribution for a given value of ^ is then obtained as a particular one-dimensional 
reduction of this two-variable distribution. The advantage of this approach is that the resulting 



skewed parton distribution satisfy automatically polynomiality conditions (5.44). However, as 
it was shown in ref. ||3l|] the parametrization of skewed quark distributions in terms of double 
distributions is not complete. This incompleteness can be seen especially clear if one considers 
the Mellin moments of skewed quark distributions. For example, the expression for H(x, A) 
in terms of double distributions requires that its A^-th Mellin moment is a polynomial of order 



A'^ — 2 in variable 4, what is in contradiction with ( p. 44 ) for even A^. This problem can be 



easily cured if one adds one additional function to the double distribution parametrization of 
light-cone nucleon matrix elements, see 

As we saw in our model the even moments of H{x,^, A^) are almost ^-independent because 
the contribution of the Dirac sea drops out. Owing to this feature of H{x, ^, A^) the additional 
function which one needs to add to the double distribution parametrization of light-cone nucleon 
matrix elements is very small. 



6 Conclusions 

We have shown that the helicity skewed distribution E is dominated in large range of A^ and ^ 



by contribution of the pion pole ( 4.36 ). This result can be viewed as generalization of well known 
chiral Ward identities for local currents to bilocal quark operators on the light cone. The fact 
that E is fixed to great extent by the pion pole contribution opens the possibility to measure H 
by choosing observables which are proportional to the product H ■ E. One of examples of such 
quantity is azimuthal spin asymmetry in hard exclusive production of pions and kaons |0] . 

The skewed quark distribution H has been computed in wide range of A^ and using 
chiral quark-soliton model. We have demonstrated that contribution of the Dirac continuum 
is crucial to describe transition between two regions \x\ > ^ and \x\ < ^. Also we saw that 
the ^ dependence of the SPD's is mostly due to Dirac continuum contribution. Since the Dirac 
continuum contribution to H^^^ is symmetric in variable x, we can expect that the dependence 

^The notion of virtual hadron can not be defined in QCD apart from special cases (large Nc limit, pions, etc.). 
We use the term "resonance part" to denote specific contributions to SPD arising only in non-forward limit. 
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of the combination H^^\x,^) — H^^\—x,S,) on ^ is rather weak. This, in particular, impUes 
that in modelhng of SPD's in terms of double distributions the function h(x, y)0 should strongly 
depend on C parity of the SPD. 

Studying the dependence of the helicity skewed parton distributions we have seen that 
the factorization ansatz H{x,(,, A'^) = ^)Ga(A^) is in contradiction with our calculations. 

For models which use the usual quark distributions to model SPD's (see e.g. ^) one 

can use replacement of the slope of Regge trajectory in small x parametrization of Aq{x) ~ 1/ 
at low normalization point by Regge trajectory oq ^ ao + Oi' Such replacement describes 
qualitatively correct the dependence of the Mellin moments of SPD's as was discussed in 
present paper. 
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A Bound-state level contribution to H^^\x, A2) and E^^\x, ^, A^) 

We present here the contributions of the discrete bound-state level to the isovector H^^\x, A^) 
and E^^^x,^, A'^). The bound-state level occurs in the grand spin K = and parity 11 = -|- 
sector of the Dirac Hamiltonian (3.11). In that sector the eigenvalue equation takes the form: 



d 



McosP(r) -|.-^-MsinP(r) | ( ho{r) 



i^,e. I 1 . (A.1) 



-M cos P(r) / \ ji(r) / \ Ji[r) 

We assume that the radial wave functions are normalized by the condition 



MsinP(r) 

or 



dr r 



hl{r)+jf{r) 



1. 



We introduce the Fourier transforms of the radial wave functions, 



(A.2) 



h{k) 



drr^ ho{r)Rkoir), 



m 



drr'^ ji{r)Rki{r), 



where 



Rkl{r) 



l^j = I [2r^ n d y sinfcr 

r '+2 ^ -K U \r dr ) r 



(A.3) 



(A.4) 



The bound-state level contribution to the H and E distribution function can be simply 
obtained from the general eqs. ( [4.19| , 4.20| ). Compactly the corresponding expressions can be 
written as: 



x\h{k)h{k')6^^ 
where 



^' A2 " 3 J (27r)2 kk' 
h{k)j{k')n^ + h{k')j{k)n'i \ + j{k')j{k) In^l + nln^ - 6^''{n^ ■ n) 



k 



k' = ^(Icx + At) +[{x-OMn-Ei, 
Unit vectors n and have components: 



n 



^(kx,(x + OMjv-Siev) 



(A.5) 

(A.6) 
(A.7) 



The individual expressions for H and E can be obtained from eq. ( [A.5 ) contracting index b with 
{6'' - ^) and Ai. 
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in the forward limit 



1.0 



1 1 

" a'=o 


1 1 1 1 1 

— Dirac continuum 
--- discrete level 
/ \ total distribution 




A /^"\ 

y / \ 




1 1 1 1 1 1 1 



-1.0 0.0 1.0 

I 

Figure 1: The isovector distribution //(x,i^,A^) in the forward hmit, A = 0. Dashed lint: 
contribution from the discrete level. Dashed- dotted line: contribution from the Dirac continuum 
according to the interpolation formula, eq. ( [4.23| ). Solid line: total distribution (sum of the 
dashed and dashed-dotted curves). 
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H''\x, \, A') 





Figure 4: Comparison of pion pole contribution and non-pole part of isovector E at various 
values of and ^. The positive curves correspond to pion pole contributions. 
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